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is strong and the flavor SU(3) singlet-octet mixing angle 6 is about zero. The effects 
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Abstract 

> ' n 

\q \ We study the n —* 77 and 77 — > -zr 77 decays using an extended three-flavor 

I Nambu-Jona-Lasinio model that includes the 't Hooft instanton induced interaction. 

We find that the //-meson mass, the n — ► 77 decay width and the n — ► 7r°77 decay 
width are in good agreement with the experimental values when the Ua(1) breaking 

OS ; 

Qh! of the Ua{X) breaking on the baryon number one and two systems are also studied. 

hi, 

*f; : 1 Introduction 

> ■ 

It is well known that the QCD action has an approximate Ul(3) x Ur(3) chiral sym- 
metry and its subsymmetry, Ua\1) symmetry, is explicitly broken by the anomaly. The 
Ua(1) symmetry breaking is manifested in the heavy mass of the 7/ meson. 

The physics of the 77 and rj' mesons have been extensively studied in the 1/Nc expan- 
sion approach [|IJ. In the Nq —* 00 limit, the £/a(1) anomaly is turned off and then the t] 
meson becomes degenerate with the pion and the 7/ meson becomes a pure ss state with 
mL{Nc — > 00) = 2m 2 K — m 2 ~ (687 MeV) 2 [0. So the ?7a(1) anomaly pushes up m v by 
about 400MeV and m v ' by about 300MeV. It means that not only the 7/ meson but also 
the 77 meson is largely affected by the anomaly. 

In order to understand the role of the Ujl(1) anomaly in the low-energy QCD, it may 
be important to study the 77- meson decays as well as its mass and decay constant. Among 
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the rj-meson decays, r] — > 77 and r] — > 7r°77 decays are interesting. They have no final 
state interactions and involve only neutral mesons so that the electromagnetic transitions 
are induced only by the internal (quark) structure of the mesons. The rj — > 77 decay 
is related to the Adler-Bell-Jackiw triangle anomaly || through the partial conservation 
of axialvector current hypothesis. For the 77 — > 7r°77 decay, it is known that the chiral 
perturbation theory (ChPT) gives too small prediction in the leading order and higher 
order terms are expected to be dominant. 

The purpose of this paper is to study the rj — > 77 and 77 — > 7r°77 decays in the 
framework of the three-flavor Nambu-Jona-Lasinio (NJL) model as a chiral effective quark 
lagrangian of the low-energy QCD. The three-flavor NJL model which involves the Z7l(3) x 
Ur(3) symmetric four-quark interaction and the six-quark flavor-determinant interaction 
0] incorporating effects of the £^4(1) anomaly is used widely in recent years to study such 
topics as the quark condensates in vacuum, the spectrum of low-lying mesons, the flavor- 
mixing properties of the low-energy hadrons, etc. In this approach the effects of the 
explicit breaking of the chiral symmetry by the current quark mass term and the {/a(1) 
anomaly on the rj — > 77 and 77 — > 7r°77 decay amplitudes can be calculated consistently 
with those on the ?7-meson mass, rj decay constant and mixing angle within the model 
applicability. 



2 Three-flavor Nambu-Jona-Lasinio model 



We work with the following NJL model lagrangian density: 



C 



£0 + £4 + C§, 

ip («<9 M 7 M — rh) if) 



G, 



a=0 



G D { det [V><(1 - 75)^'] + det [^(1 + 7s)^] } 



(1) 
(2) 

(3) 
(4) 



Here the quark field tp is a column vector in color, flavor and Dirac spaces and X a (a = 
... 8) is the U(3) generator in flavor space. The free Dirac lagrangian Cq incorporates the 
current quark mass matrix rh = diag(m u , rrid, m s ) which breaks the chiral Ul(3) x Ur(3) 
invariance explicitly. £4 is a QCD motivated four-fermion interaction, which is chiral 
Ul(3) x Ur(3) invariant. The 't Hooft determinant £q represents the ^(1) anomaly. It 
is a 3 x 3 determinant with respect to flavor with i,j = u, d, s. 

Quark condensates and constituent quark masses are self-consistently determined by 
the gap equations in the mean field approximation. The covariant cutoff A is introduced 
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to regularize the divergent integrals. The pseudoscalar channel quark-antiquark scattering 
amplitudes are then calculated in the ladder approximation. From the pole positions of 
the scattering amplitudes, the pseudoscalar meson masses are determined. We define the 
effective meson-quark coupling constants g mq and g wqq by introducing additional vertex 
lagrangians, 

£ W s = g mq ipi^ v ^<P v ^ ( 5 ) 

£nqq = #^^75 A 3 ?/^ , (6) 

with X n = cos 9X 8 — sin ^A . Here <fi is an auxiliary meson field introduced for convenience 
and the effective meson-quark coupling constants are calculated from the residues of 
the gg-scattering amplitudes at the corresponding meson poles. Because of the SU(3) 
symmetry breaking, the flavor A 8 — A components mix with each other. Thus we solve 
the coupled-channel qq scattering problem for the 77 meson. The mixing angle 9 is obtained 
by diagonalization of the gg-scattering amplitude. It should be noted that 9 depends on 
q 2 . At q 2 = m 2 , 9 represents the mixing angle of the A 8 and A components in the 
77-meson state. In the usual effective pseudoscalar meson lagrangian approaches, the 77 
and rf mesons are analyzed using the g 2 -independent 77-77' mixing angle. Because of the 
g 2 -dependence, 9 cannot be interpreted as the 77-77' mixing angle. The origin of the q 2 - 
dependence is that the 77 and 77' mesons have the internal quark structures. The meson 
decay constant J'm (M = tt, K, 77) is determined by calculating the quark-antiquark one- 
loop graph. The explicit expressions are found in ||. 

3 77 — > 77 decay amplitude 

The 7T°, 77 — >• 77 decay amplitudes are given by 

(7(^)7(^2)^(5)) = t{2'Kf5\k 1 + k 2 - q)e^ pa e^k°f M ^(q 2 ) , (7) 

where e\ and e 2 are the polarization vectors of the photon. By calculating the pseudoscalar- 
vector-vector type quark triangle diagrams, we get the following results. 

7; _ 77 = -g n F(u, tt ) , (8) 
air 

T v ^ n = -g^—^ cos9 {5F(u,r)) -2F(s,r})} 



sm9V2{5F(u,ri) + F(s,r])}} . (9) 



Here a is the fine structure constant of QED and F(a, M) (a = u,s and M = 7r°, 77) is 
defined as 

F(a, M) = ( l dx I' dy — ^ ~^ M ° r . (10) 

v ' ; Jo Jo y Ml -m 2 M x(l -x){\ -y) v ; 
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Then the M — > 77 decay width T(M — > 77) is given by Y(M — > 77) = |7M^ 77 | 2 m^/(647r). 

In the chiral limit, the pion mass vanishes and F(u, 7r°) becomes 1/M U . In this limit, 
the Goldberger-Treiman (GT) relation at the quark level, M u = g n f n , holds in the NJL 
model and this leads to 7^.o^ 77 = a/(7r/ 7r ) which is same as the tree-level results in the 
Wess-Zumino-Witten lagrangian approach [0. It should be mentioned that we have to 
integrate out the triangle diagrams without introducing a cutoff A in order to get the 
above result though the cutoff is introduced in the gap equations in the NJL model. In 
the U(3)l x U(3)r version of the NJL model, the WZW term has been derived using the 
bosonization method with the heat-kernel expansion [§, [J. In their approach, 0(1/ A) 
term has been neglected and it is equivalent to taking the A — > 00 limit. 



4 77 —> 7r°77 decay amplitude 

The r] — > 7r°77 decay amplitude is given by 

(^(pMh, e l ) 1 (k 2l e 2 )\ V ( P )) = i(27r) 4 5\p n +h + k 2 - p)e^ 2 T, u . (11) 

The dominant contributions to this process in this model are the quark-box diagrams. 
Following the evaluation of the quark-box diagrams performed in JTDJ , we obtain 



with 



I r d q ^ 

Tyv = -i-j=(cos6 - V2 sin 9)e 2 g mq g« qq J j^gJ2 U l^ ( 12 ) 



fi-M + te l °fi + p , -k' l -k f 2 -M + te 

X lu i + i> - % - M + ie lfl i + i> - M + te } ' (13) 
u »v = tr 175i ~, , - 7s- 



- M + ie i + fa- M + ie 
X lu i + i)-l/ l ~M + le 1 ^ + $ - M + le } ' (14) 



i lo i-M + ie '"t + fc-M + ie 

1 1 1 

X ^^/i + ^ + h-M + ie^i + ^-M + it) 1 (15) 

K = U^k^h), (16) 

U% = U 2 ^(h^k 2 ), (17) 

U% = U 3 ^(h^k 2 ). (18) 
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Here M is the constituent u,d-quark mass. Because the loop integration in QI2] ) is not 
divergent, we again do not use the UV cutoff. Then the gauge invariance is preserved. 
The inclusion of the cutoff that is consistent with the gap equation will break the gauge 
invariance and make the present calculation too complicated. Note that the strange quark 
does not contribute to the loop. 

On the other hand the amplitude T MJ , has a general form required by the gauge invari- 
ance ir ' ' ' 
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= A(x 1 ,x 2 )(k^-k 1 -k 2 g 



+ B(x h x 2 ) 



k\ ■ k 2 „ v . a v „ . v 

^—pp + x\k 2 p + x 2 p^k, 

mi 



(19) 



with Xi = p - ki/rrii. With A and B, the differential decay rate with respect to the energies 
of the two photons is given by 



d 2 Y 

dx\dx 2 



mi 



256tt 2 



1 2 
A + -B 
2 



mz 



1 2 



2( Xl +x 2 ) + ^-l 



mi 



+ ]\ B \ 2 
4 



4xiX 2 



mz 



2( Xl + x 2 ) + ^-l 



nif 



(20) 



Though the mass of r\ as a qq bound state depends on G e p, we use the experimental value 



m. 



and 



547 MeV in evaluating 

1 

2 



). The Dalitz boundary is given by two conditions: 

,2 \ 



1 - 



mz 



m 



< X\ + x 2 < 1 - 



m 



v 



xi + x 2 - 2x x x 2 < - I 1 



mz 



(21) 



(22) 



In evaluating ([13D - ([T8|) , one only has to identify the coefficients of p^p v and g^ v . Details 
of the calculation are given in [JTU]. Defining A and B by 



d A a 6 



(2vr) 4 



-% Ag^ + B 



i=l 



ml 



we find A and B as 
A = 



B 



1 2 
-^=(cos6» - V2sin6)e 2 g 7Tqq g mq 

—=(cos9 — \Pl si 
a/3 



m 2 cr L 
2 B 

smv)e~g mq g mq — -, 



B 

A — 2xix 2 — 
a 



with 



a 



(h + k 2 f 



mz 



2( Xl + x 2 ) + -%-l. 



m 2 m 2 



(23) 

(24) 
(25) 

(26) 



We evaluate A and B numerically and further integrate (^) to obtain the r\ — > 7r°77 decay 
rate. 
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5 Numerical Results 



The recent experimental results of the tt°,i] — > 77 decay widths are r(7r° — > 77) 



7.7 ± 0.6 eV and T(r] — > 77) = 0.510 ± 0.026 keV |L2| and the reduced amplitudes are 



+ 77 



T 

r \~ *77 



(2.5 ±0.1) x 10~ n [eV] -1 , (27) 
(2.5 ± 0.06) x 10~ n [eV]" 1 . (28) 



From Eq. @ and Eq. (Q), we get 7^_^ 77 = (5/3)7^o^ 77 in the U "^(1) limit. Therefore 
in order to reproduce the experimental value of 7^_^ 77 , the effect of the ?7a(1) anomaly 
should reduce 7^_^ 77 by a factor 3/5. On the other hand, the experimental value of the 
i] — » 7r°77 decay width is 



1^(77^7^77) = 0.85 ±0.19 eV. (29) 

In our theoretical calculations, the parameters of the NJL model are the current quark 
masses m u = m^, m s , the four-quark coupling constant G$, the six-quark determinant 
coupling constant Gd and the covariant cutoff A. We take as a free parameter and 
study 7] meson properties as functions of Gd. We use the light current quark masses 
mu = m d = 8.0 MeV to reproduce M u = M d ~ 330 MeV (~ 1/?>M N ) which is the value 
usually used in the nonrelativistic quark model. Other parameters, m s , Gd, and A, are 
determined so as to reproduce the isospin averaged observed masses, m n ,mK, and the 
pion decay constant f^. 

We obtain m s = 193 MeV, A = 783 MeV, M u4 = 325 MeV and g nqq = 3.44, which 
are almost independent of Gd- The ratio of the current s-quark mass to the current u,d- 
quark mass is m s /m u = 24.1, which agrees well with m s /m = 25±2.5 (m = |(m u ±md)) 
derived from ChPT The kaon decay constant fx is the prediction and is almost 

independent of Gd- We have obtained fx = 97 MeV which is about 14% smaller than 
the observed value. We consider this is the typical predictive power of the NJL model in 
the strangeness sector. 

Table 1 summarizes the fitted results of the model parameters and the quantities 

necessary for calculating the rj — > 77 and r] — > 7r°77 decay widths which depend on Gd- 

We define dimensionless parameters G<§ = -G d (A/2tt) 4 AA^ c 2 and Gf = G s (A/2n) 2 N c . 

When G e p is zero, our lagrangian does not cause the flavor mixing and therefore the ideal 

mixing is achieved. The "77" is purely uu + dd and is degenerate to the pion in this limit. 

We next discuss the tt° — > 77 decay. The calculated result is 7^o^ 77 = 2.50 x 
-11 



10 (1/eV) which agrees well with the observed value given in Eq. (|27|) . The current 
algebra result is T^o^,^ = a/(irf n ) = 2.514 x 10~ n (l/eV), and thus the soft pion limit is 
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Table 1: The parameters of the model, the rj — ► 77 decay amplitude 7^_> 77 and 77 
7r°77 decay width r(77 — > 7r°77) for each Gff 





Gf 


M s [MeV] 


m v [MeV] 


9 [deg 




9qqq 


T 

V~ >77 


1/eV] 


T(rj -> vr°77) [eV] 


0.00 


0.73 


556 


138.1 


-54.74 


3.44 


4.17x 


10- 11 


2.88 


0.10 


0.70 


552 


285.3 


-44.61 


3.23 


3.95x 


10- 11 


2.46 


0.20 


0.66 


549 


366.1 


-33.52 


3.12 


3.68x 


10- 11 


2.06 


0.30 


0.63 


545 


419.1 


-23.24 


3.11 


3.39x 


10- 11 


1.71 


0.40 


0.60 


541 


455.0 


-14.98 


3.15 


3.10x 


10- 11 


1.42 


0.50 


0.57 


537 


479.7 


-8.86 


3.20 


2.86x 


10- 11 


1.20 


0.60 


0.54 


533 


497.3 


AAA 

-4.44 


3.25 


2.65x 


10- 11 


1.04 


0.70 


0.51 


529 


510.0 


-1.25 


3.28 


2.48x 


10- 11 


0.92 


0.80 


0.47 


525 


519.6 


1.09 


3.30 


2.35x 


10- 11 


0.84 


0.90 


0.44 


522 


527.0 


2.84 


3.31 


2.23x 


10- 11 


0.77 


1.00 


0.41 


518 


532.8 


4.17 


3.32 


2.14x 


10- 11 


0.71 


1.10 


0.40 


514 


537.5 


5.21 


3.32 


2.07x 


10- 11 


0.67 


1.20 


0.35 


511 


541.3 


6.02 


3.31 


2.00x 


10- 11 


0.63 


1.30 


0.32 


507 


544.5 


6.66 


3.30 


1.95x 


10- 11 


0.61 


1.40 


0.29 


504 


547.2 


7.17 


3.29 


1.91x 


10- 11 


0.58 


1.50 


0.25 


500 


549.4 


7.57 


3.28 


1.86x 


10- 11 


0.56 


1.60 


0.22 


497 


551.4 


7.90 


3.26 


1.83x 


10- 11 


0.55 



a good approximation for 7r° — > 77 decay. The chiral symmetry breaking affects 7^ r o_ n , 7 in 
two ways. One is the deviation from the G-T relation and another is the matrix element 
of the triangle diagram F(u, 7T°). Our numerical results are g n = 3.44, M u j f n = 3.52 and 
F(u, 7r°)M u = 1.015, therefore the deviations from the soft pion limit are very small both 
in the G-T relation and the matrix element of the triangle diagram. 

Let us now turn to the discussion of the rj — > 77 decay. The calculated results of 
the f] — > 77 decay amplitude 7^ 77 are given in Table 1 and shown in Fig. 1. The 
experimental value of the rj — > 77 decay amplitude is reproduced at about G e ^ = 0.7. 
The calculated ?7-meson mass at Gff = 0.7 is m v = 510 MeV which is 7% smaller than 
the observed mass. G e § = 0.7 corresponds to Gd{ss) /Gs = 0.44, suggesting that the 
contribution from £ 6 to the dynamical mass of the up and down quarks is 44% of that 
from £4. 

The mixing angle at G*§ = 0.7 is 9 = —1.3° and that indicates a strong OZI violation 
and a large (u,d)-s mixing. This disagrees with the "standard" value 9 ~ —20° obtained 
in ChPT[14]. This is due to the stronger {7a(1) breaking in the present calculation. The 
difference mainly comes from the fact that the mixing angle in the NJL model depends on 
q 2 of the qq state and thus reflects the internal structure of the 77 meson. On the contrary 
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Figure 1: Dependence of the 77 — > 77 decay amplitude on the dimension-less coupling 
constant G c § . The horizontal dashed line indicates the experimental value. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



D 



the analyses of ChPT[14] assume an energy-independent mixing angle, i.e., 8(q = m 2 ) = 
9(q 2 = m 2 ,). 

The 77 decay constant is almost independent of Gd- f v = 91.2 MeV (~ /„-) at G c § = 
0.7. Therefore it seems that the 77 meson does not lose the Nambu-Goldstone boson nature 
though its mass and mixing angle are strongly affected by the £7a(1) breaking interaction. 

The 77 — > 7r°77 decay gives us an independent information on the structure of the 7? 
meson. The calculated 77 — > 7r°77 decay widths are also given in Table 1 and shown in 
Fig. 2. At G c p = 0.70, we obtain T(r] — > 7r°77) = 0.92, which is in good agreement 
with the experimental data shown in (|2D|). This process was studied in ChPT|T3[ and 
in the extended NJL model |TB[. The difference between these approaches and ours are 
discussed in 



6 Effects of the £/a(1) Anomaly in Baryons 

Since the effects of the Ua(1) anomaly are rather large in the pseudoscalar sector, 
it is natural to ask if one can see some effects in the baryon sector. It was pointed 
that the instanton can play an important role in the description of spin-spin 



out m 

forces, particularly for light baryons. The pattern of these effects can be very hard to 
disentangle from one-gluon exchange. The effects of the instanton induced interaction 
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Figure 2: Dependence of the r] — > 7r°77 decay width on the dimension-less coupling con- 
stant G*§ . The horizontal solid line indicates the experimental value and the dashed lines 
indicate its error widths. 

T(rj tt°77) [eV] 

3 i i i i i i i i i i i i i i i i i 




in baryon number B = 2 systems were studied in |L9|. It was shown that an attraction 
between two nucleons is obtained by the two-body instanton induced interaction, while 
the three-body interaction is strongly repulsive in the H-dibaryon channel and makes the 
H-dibaryon almost unbound. 

We estimate the effects of the £/a(1) anomaly on the B = 1 and B = 2 systems by 
employing the six-quark determinant interaction given in Eq. (D) whose strength was de- 
termined so as to reproduce the observed ^-meson mass, the r] — > 77 decay width and the 
7] — > 7r°77 decay width, namely, Gff = 0.7. It is done by calculating the matrix elements 
of the the Ua{^-) breaking interaction hamiltonian with respect to unperturbed states of 
the MIT bag model and the nonrelativistic quark model (NRQM). For B = 2 systems, 
we only consider the (OS) 6 configuration of the six valence quark states. Therefore, the 
matrix element with respect to such a state gives a measure of the contribution of the 
[/^(l) breaking interaction either to the dibaryon or to the short-range part of the inter- 
action between two baryons. The determinant interaction induces not only three-body 
but also two-body interactions of valence quarks when the vacuum has a nonvanishing 
quark condensate. The details of the calculation are described in [ pOf . 

Table 2 shows the contribution of the two-body term for B — 1. The contribution 
to the decuplet baryons vanishes in the SU(3) limit and therefore comes only from the 
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Table 2: Contribution of the two-body term to octet and decuplet baryons. All the entries 
are in units of MeV 



wave function 


N E E A 


AS* E* n 


MIT 

NRQM 


-43.9 -41.2. -41.2 -42.9 
-40.88 -36.6 -36.6 -39.4 


0.12 0.12 
0.07 0.07 



Table 3: Baryon component, SU(3) multiplet, spin, isospin and strangeness of the eight 
channels of two octet baryons 



channel 


Baryon component 


SU(3) multiplet 


Spin 


Isospin 


Strangeness 


I 


NN 


10* 


1 








II 


NN 


27 





1 





III 


NT 


27 





3/2 


-1 


IV 


NT-NA 


27 





1/2 


-1 


V 


iVE — iVA 


10* 


1 


1/2 


-1 


VI 


NT 


10 


1 


3/2 


-1 


VII 


NT — NA 


8 


1 


1/2 


-1 


VIII 


H 


1 








-2 



S77(3) asymmetry of the quark wave function. The three-body term does not contribute 
to the B = 1 states. Thus the NA mass difference due to the ?7a(1) breaking interaction 
is about 15% of the observed one. 

We next discuss the case of B = 2. We consider all the possible channels which 
are made of two octet baryons listed in Table 3. Table 4 shows the contribution of the 
two-body term. The channel VIII gets the strongest attraction, about 170 MeV, and the 
channel VII gets the second strongest attraction. The contributions of the three-body term 
to the H-dibaryon and strangeness —1 channels are given in Table 5. It should be noted 
that the three-body term has no effect on the NN channels, and that the contributions to 
the channels III, IV and V reflect the SU(3) breaking in the quark wave function. The 
contributions of the three-body term in channels VI, VII and VIII are remarkable and 
one will be able to observe some effects experimentally. 

We should comment on the difference between the determinant interaction used here 
and the instanton- induced interaction used in ref. [T9]. The relative contributions of the 
Ua(1) breaking interaction within the baryonic sector or within the mesonic sector are 
similar for the two interactions. However, the ratio of those in the baryonic sector to 
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Table 4: Contributions of the two-body term to the eight channels of two octet baryons 
listed in Table 3. All the entries are in units of MeV 



wave function 


I 


II 


III 


IV 


V 


VI 


VII 


VIII 


MIT 


-89.9 


-85.1 


-88.9 


-86.3 


-93.9 


-96.5 


-121.5 


-162.6 


NRQM 


-120.2 


-105.3 


-102.3 


-103.9 


-118.3 


-116.7 


-148.0 


-182.6 



Table 5: Contributions of the three-body term to the H-dibaryon (VIII) and strangeness 
— 1 two octet baryon channels (III- VII). All the entries are in units of MeV 



wave function 




III 




IV 




V 


VI 


VII 


VIII 


MIT 


-6 


x lO" 2 


-6 


x 10~ 2 


-7 


X 10" 2 


20.7 


25.1 


40.7 


NRQM 


-5 


X ID" 2 


-5 


X ID" 2 


-5 


X ID" 2 


28.3 


34.9 


56.1 



those in the mesonic sector is about |. Namely, if one fixes the strength of the interaction 
so as to give the same mass difference of 77 and rj', the effects of the instanton-induced 
interaction in the baryonic sector would be about | stronger than those of the determinant 
interaction. After this correction the strength of the present Ua{1) breaking interaction is 



consistent with that used in the calculation of the baryon-baryon interaction in ref. [19 



7 Summary 

We have studied the 77 — > 7r°77 decay in the three-flavor NJL model that includes the 
[/^(l) breaking six-quark determinant interaction. The 77 meson mass, the 77 — > 77 decay 
width and the 77 — > 7T°77 decay width are reproduced well with a rather strong Ua{^-) 
breaking interaction, that makes rji — rj s mixing angle 6 ~ 0°. The effects of the Ua{^-) 
breaking on the baryon number one and two systems have been also studied. 

Finally, we should note that the NJL model does not confine quarks. Since the NG 
bosons, 7r, K, and 77, are strongly bound, the NJL model can describe their properties 
fairly well. However the 7/ meson state in the NJL model has an unphysical decay of the 
77' — > qq. Therefore we do not apply our model to the 7/ meson. Further study of the 
[7^4 (1) breaking will require the construction of a framework which can be applied to the 
77' meson. 
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